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Abstract. We study the phase structure of a surface model by using the canonical Monte Carlo simulation 
technique on triangulated, fixed connectivity, and spherical surfaces with many fine holes. The size of 
a hole is assumed to be of the order of lattice spacing (or bond length) and hence can be negligible 
compared to the surface size in the thermodynamic limit. We observe in the numerical data that the model 
undergoes a first-order collapsing transition between the smooth phase and the coUapsed phase. Moreover 
the Hasudorff dimension H remains in the physical bound, i.e., H < 3 not only in the smooth phase but also 
in the collapsed phase at the transition point. The second observation is that the collapsing transition is 
accompanied by a continuous transition of surface fluctuations. This second result distinguishes the model 
in this paper and the previous one with many holes, whose size is of the order of the surface size, because 
the previous surface model with large-sized holes has only the collapsing transition and no transition of 
surface fluctuations. 
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1 Introduction 



Biological membranes such as the so-called cell membranes 
or plasma membranes have many holes called transport 
protein or protein channel, through which some biologi- 
cal materials are transported in/out. Currently it is well 
known that cell membranes are heterogeneous due to cy- 
toskeletons and membrane proteins including the trans- 
port proteins, and moreover, the membranes have complex 
structures such as rafts and fences [1,2 . 

The fluid mosaic model [3] seems to be valid only when 
these structures of size smaller than the membrane size 
were neglected. By neglecting further the structures of 
size negligible compared to the membrane size, we have 
a homogeneous surface, which has been described by the 
conventional surface model of Helfrich and Polyakov [31[51 
[6] . The crumpling phenomena was extensively studied the- 
oretically and numerically [7u8^i9l i10u11|. i12 ^ i13| . i14J by using 
this homogeneous surface model. Current understanding 
of membranes on the basis of statistical mechanics are re- 
viewed in p!51[TMT7[fT5] . We also know that the collapsing 
transition is accompanied by a transition of surface fluc- 
tuations, and both transitions are of first-order [l9 l l20 l [2T ] . 

It was recently reported that a surface model with 
many holes undergoes a collapsing transition between the 
smooth phase and the collapsed phase, and no transition 
of surface fluctuations can be seen in that model [52] . The 
results in |22j are considered to be reflecting an inhomoge- 
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neous structure in biological membranes. The Haniiltonian 
of the model in [22 is the one of Helfrich and Polyakov. 



However, the size of holes in [22] is assumed to be of 
the order of the surface size in the limit of A^ ^ oo, {N is 
the total number of vertices) or in other words in the ther- 
modynamic limit. Thus, the phase structure of the surface 
model with small sized holes still remains to be studied. 
We expect that the phase structure of the homogeneous 
model is influenced by such holes, because the vertices at 
the edge of the holes are relatively freely moving; no bend- 
ing energy is defined on the edges, i.e. those vertices are 
considered to be in a free boundary condition. For this 
reason, no transition of surface fluctuations can be seen 
in the model in [22]. Thus, the influence of holes on the 
phase structure is expected to persist in such a case that 
the size of holes reduces sufficiently small compared to the 
surface size. 



In this paper, we study a surface model of Helfrich and 
Polyakov on triangulated, fixed connectivity and spherical 
surfaces with many fine holes by Monte Carlo (MC) simu- 
lations. The purpose of the present paper is to see whether 
or not the phase transitions of the conventional homoge- 
neous model are influenced by the presence of many fine 
holes and change their properties. 
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2 Model 



The partition function we are studying is of the form 



./ N 

Z= Y[dX,exp[-S{X,T)], 



(1) 



4=1 



SiX,T)^Si + bS2. 



The Hamiltonian S is defined on triangulated lattices. The 
construction technique for the lattice with holes will be 
presented below. The symbols X and T in S* denote the 
vertex position and the triangulation, where T is fixed in 
the simulations. / denotes that the center of mass of the 
surface is fixed. The Gaussian bond potential Si and the 
bending energy S2 are defined so that 



Si = J2ix, - x,f 

(ij) 



s. 






n,), (2) 



where the variable Xi{^ R'^) denotes the position of the 
vertex i, and ni(G S^= the unit sphere in R'^) denotes a 
unit normal vector of the triangle i. 

The construction technique for the triangulated lat- 
tice is almost identical to that for the lattices in [22] and 
is summarized as follows: We start with the icosahedron, 
which is characterized by 12 vertices, 30(= 3x12—6) bonds, 
and 20(=2xl2— 4) triangles. By dividing the bonds of the 
icosahedron into d pieces of uniform length a, we firstly 
have a triangulated spherical surface of size iVo = 10£^-|-2, 
which is the total number of vertices on the surface with- 
out holes. Secondly, a sublattice is obtained by dividing t 
edges into m partitions (m = 1, 2, ■ • •), where each parti- 
tion is of leneth L = f.lm, in the unit of a if m, divides I. 

□□ 




Fig. 1. Surfaces with holes of size (a) (iV, L) = (2942, 3) and 
(b) (A'^, L) — (3402, 4), which can also be characterized by two 
integers {£,m) = (18,6) and (£, m) = (20, 5), respectively. The 
holes in (a), (b) respectively contain 6, 13 triangles, which are 
excluded from the surface. 



We should note that the size of holes can be charac- 
terized by L. Two types of lattices corresponding to 



L = 3, L = 4, 



(3) 



are assumed in this paper. The minimum size of holes is 
given by the condition L = 3, where the center of hole 
is a vertex that is excluded from the surface. The sec- 
ond minimum size is given by L = 4, where each hole 
has 3 vertices that are excluded from the surface. Figures 
IHa) and [2b) show the lattices of {N,L) = (2942,3) and 
{N, L) = (3402, 4), which are given by {(., m) = (18, 6) and 
{£,m) = (20,5), respectively. The size of holes given by 
Eq. ([3]) remains fixed while the total number of vertices iV 
is changed, and therefore the size of holes becomes negli- 
gible compared to the surface size in the thermodynamic 
limit in both cases in Eq.((3]). 

The reason why we assume these two values in Eq. (|3|) 
for L is see that the final results obtained from the model 
are independent of L. In fact, both of the assumed L are 
negligible to the surface size in the thermodynamic limit, 
and hence the results should not be independent of L in 
Eq.([3l). 

The total number of holes in one face of the icosahe- 
dron is given by m{m — l)/2 and, hence the total number 
of holes over the surface is 10TO(m— 1) . Because of the holes 
on the surface, the total number of vertices N of the lattice 
with holes are reduced from A^'o- (i— 2)(L— 1)/2(=^^ ~ i) 
vertices per a hole are excluded from the surface. Then, 
we have N — No~5{m — l){£ — in){£~2m)/m, because N 
is given by No — lOm{m — l) x (L — 2)(L — 1)/2, which can 
also be expressed as N — No~5{m—l){£ — m){£ — 2m)/m 
by using L = i/m. The lattice size is given by two inte- 
gers {£, m), and hence the sizes N and L are expressed by 
(iV, L) = {10i'^+2-5{m-l){i-m){£-2m)/m, £/m).We note 
that N includes the vertices on the boundary of holes. 

The ratio Rl{oo) of the area of holes to that of the sur- 
face including the holes is given Rl{oo) -^ (i^ — 3)/{2L^) 
in the limit of tti— >oo and ^— >cx3 under constant L~£/m. 
In fact, the total number of triangles in a hole is L^ — 3 
and then, the total number of triangles in the holes is 
10TO(m — 1)(L^ — 3), which is easily understood since the 
total number of faces in the icosahedron is 20, and the 
total number of holes in a face is m{m~\)/2 as stated 
above. On the other hand, the total number of trian- 
gles on the triangulated sphere is 2iVo — 4. Then, we have 
Rl{N) = 10m(m-l)(L2_3)/(2iV(r4). By using L^ = {£/m)'^ 
and £^ = (iVo — 2)/10 and identifying the factor {m—l)/m 
in Rl{oo) as 1, we have Rl{oo) = (L^ — 3)/(2L^) in the 
limit of both ?7i— >cxd and £— »oo with finite L = £/m. 

We show in Table [T] some of the numbers that char- 
acterize the lattices we use in the MC simulations. In the 
case of the model in [35], in is fixed and then the ratio 
Rm{oo) of the area of holes to that of the surface includ- 
ing the holes is given i?m(oo) = (rn — l)/(2m) in the limit 
of A^'o -^ cxo. In that case in [55] Rm{N) is almost identical 
with Rra{oo) even when N is relatively small A^(~ 10^), 
while Rl{N) in this paper deviates from Rl{oo) about 
10% even on the largest surfaces. The reason for the devi- 
ation of Rl{N) from Rl{oo) is the constraint m—KX). In 
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Table 1. The total number of vertices N and N 
m, and the corresponding Rm{N). No is the total 
vertices of the original lattice, which has no holes. 



i£,m) 


(iV,L) 


No 


-Rl(oo) 


(18,6) 


(2942,3) 


3242 


l/3(~ 0.333) 


(24,8) 


(5202,3) 


5762 


1/3 


(30,10) 


(8102,3) 


9002 


1/3 


(39,13) 


(13652,3) 


15212 


1/3 


(20,5) 


(3402,4) 


4002 


13/32(~ 0.406) 


(28,7) 


(6582,4) 


7842 


13/32 


(36,9) 


(10802,4) 


12962 


13/32 


(44,11) 


(16062,4) 


19362 


13/32 



fact, the number m is relatively small compared 
expression Rl{co) deviates from Rl{N) just by 
For this reason, the finite size effect influences i 
in this paper more strongly rather than the moc 

We note that Rm{oo) in [22] is 1/3 and 3/E 
while Rl{oo) in this paper is 1/3 and 13/32. [ 
the two ratios assumed in this paper are both ah 
tical to those in [22J . The only difference betwee 
models in this paper and in [22' is in the size oi 

The canonical Metropolis technique is emplo; 
date the variable X of the model. The variable 
domly shifted to a new position X' such that X' 
where SX is a position in a small sphere with ra 
at the beginning of the simulations to maintain a 
acceptance rate. The new position X' is acceptec 
probability Min[l,exp(— (55)], where SS — S{ne^ 



□ 




Fig. 2. Snapshots of the surface {N,L) = (13652,3) in (a) 
the collapsed phase and (b) the smooth phase at the transition 
point b— 1.43. (c) The surface section of (a), and (d) the surface 
section of (b). The mean square size X^ is about X'^ — 63 in 
(a), and X^^ 141 in (b). 



3 Results 

First we show snapshots of surfaces in Figs.[2la) Q 
and the corresponding surface sections in Figs. Q 
[Sfd). These were drawn in the same scale and ouiainea 
at 6 = 1.43, which is the transition point of the surface 
{N, L) = (13652, 3). We should note that two phases, the 
collapsed and the smooth phases, coexist even on such 
large surface at &=1.43. The snapshots in Figs. [SI a) and 
[DJb) are typical of the collapsed phase and the smooth 
phase at the transition point b — 1.43. The mean square 
size X'^ of the surfaces in Figs.[2l[a) and[2Ub) is given by 
X^ = 63 and X"^ = 141, respectively, where X'^ is defined 
as follows: 



X' 






X] 



X 



N ^ 



(4) 



where X is the center of mass of the surface. 

The mean square size X'^ versus b are shown in Figs. 
[31Ja) and[3]Jb). We see that the transition point be moves 
right with increasing N just like in the model of [23- The 
value of be is relatively larger than that &c — 0.77 of the 
model without holes in [20^ , and b^. of the model in this 
paper is also almost identical to that of the model with 
many large-sized holes in }22j . 
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Fig. 3. The mean square size X versus h of the surfaces (a) 
L = 3 and (b) L = 4. 
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The variance (7^2 oi X^ is defined by 

C^.^^{{X^-{X^)f)^ 



(5) 



and then size fluctuations are expected to be reflected 
in Cx2 . Figures HJ^a) and SJ^b) show Cx'^ versus b. The 
anomalous peaks seen in Cx^ imply large fluctuations of 
the surface size and indicate a collapsing transition, just 
like in the cases of the model without holes f20l and the 
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0.08 




0.04 



Fig. 4. The variance Cx^ 
and (b) L = 4. 



versus b of the surfafp« fa^ T, 



It is interesting to see whether the collapsed phase is 
physical in the sense that the collapsed surface can appear 
in the three-dimensional space. We know that collapsed 
surfaces of the model with/ without holes are characterized 
by Hausdorff dimension _ff < 3 at the collapsing transition 
point, i.e., the collapsing transition is physical even though 
the models are just the so-called phantom surface model 
because of the self-intersecting property. To see the Haus- 
dorff dimension in the collapsed phase at the transition 
point be, we firstly show in Figs. [6ja)-[6l^f ) the variation of 
X^ against MCS (Monte Carlo sweeps) at be of the largest 
three surfaces in both cases L — i and L — A. Four dashed 

and 



lines in each figure denote the lower bound X^ 

, -.^2 r.oUsmn') ^ 



2 col(sinoJ 
mill 



model with large-sized holes [22] . The wori 
denotes the property that the peak value go 
C-max ^ oo in the limit of iV ^ oo. This ] 
actually be confirmed as a scaling property. 



-X2 




0.05 



5000 10000 



(a) 



N 



(b) 



5000 



™^'' no- 

x2 a,g< 



Fig. 5. Log-log plots of the peak value C 

cases (a) L = 3 and (b) L = 4. The straight lines -v 

fitting the data to Eq.(l6|. 



The peak values C'^2^ are plotted in a 
against A'' in Figs.[5ja) andtS^b). The scalinj 
C™2^ can be seen in the fit of the form 



^X2 



N", 



where z/ is a critical exponent of the coUapsir 
The straight lines in the figures are obtained 
data to Eq.® with the results 

i'L=3 = 0.97 ± 0.12 (L = 3), 
1/^^4 = 1.01 ±0.07 (L = 4). 

The finite-size scaling theory indicates that t 
is of first-order in both of the cases L = J 
because both of the exponents in Eq.([7]) ar 
The fact that i'l=3 and 1/^=4 are almost the 
the errors is consistent to the expectation t 
be independent of L. 



□ 



4 xl(f 

sweeps 




4 S xlO 

sweeps 



4 8 slQp 
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Fig. 6. The variation of X^ against MCS at (a)-(c) L = 3, 
and (d)-(f) L = 4. The distribution h{X'^) (histogram) of X^ 
on the surfaces (g) L^3,N = 5202 and (g) L = 4,N = 6582. 
The data were obtained at the transition point bc{N), which 
depends on both N and L. 



In Figs.JHKg) andini^h), we plot the distribution h{X'^) 
(histogram) of X^ obtained at the transition points of the 
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surfaces L = 3,N = 5202 and L = 4:,N^ 6582. Both of the 
histograms have a double peak structure, which clearly 
indicates a first-order collapsing transition. The double 
peak structure can also be expected on larger surfaces 
such as A'' = 13652 and N = 16062, however, it is more 
hard to see the double peak in the histogram on such 
large surfaces. 



Table 2. The lower bound X^f^' and the upper bound X^^°^ 
in the collapsed state, and the lower bound X^™° and the 
upper bound X^™° 
the smooth state. 



for obtaining the mean value X (smo) in 



iV 



x„ 



x„ 



x„ 



x„ 



3 


2942 


1.34 


12 


24 


27 


44 


3 


5202 


1.39 


20 


44 


48 


76 


3 


8102 


1.41 


27 


55 


62 


108 


3 


13652 


1.42 


35 


93 


100 


165 


4 


3402 


1.44 


15 


35 


39 


59 


4 


6582 


1.51 


28 


58 


65 


106 


4 


10802 


1.56 


29 


95 


105 


170 


4 


16062 


1.58 


60 


140 


150 


232 



where H is the Hausdorff dimension. We have i/smo of the 
smooth phase and Hco\ of the collapsed phase such that 

iJsmo = 2.30 ± 0.27, iJcoi = 2.62 ± 0.60 (L = 3), 
iJsmo == 2.16 ± 0.19, iJcoi = 2.30 ±0.41 (L = 4). (9) 

We find that i/smo in both cases L = 3 and L = 4 is almost 
identical to the topological dimension H = 2 as expected, 
and moreover that the collapsed phase is also considered 
to be physical because ifcoi is less than 3 although the 
errors are relatively large. 
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The lower bound X^^^ ^ ' and the upper bound X, 
and some other numbers characterizing the transition are 
shown in Table [21 The symbol he in Table [2] denotes the 
bending rigidity where the mean values X^ are computed 
by using the upper and the lower bounds. 
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Fig. 7. Log-log plots of the mean values X"^ against A^, where 
X^ are obtained by averaging the series of X^ with X^i„, X^^^ 
shown in Table [2l 



In Figs. Efa) and Efb) we show the mean values X^ 
against A^ in a log-log scale, where X^ are obtained by 
using the variations of X^ and the upper and the lower 
bounds at he shown in Table [2J The error bars in Figs. 
[Tlja) andlljb) denote the standard deviations. The straight 
lines drawn in the figure are obtained by fitting the data 
to the scaling relation 



(smo) 

Fig. 8. X^ obtained by undercooling and superheating the 
surfaces of (a) L = 3,iV = 8102 and (b) L = 4, TV = 6582. 



Figures[5I^a) and[5I^b) show a hysteresis of X^ obtained 
in the undercooling and superheating processes on the sur- 
faces of L = 3, iV = 8102 and L = 4, TV = 6582, where 
the undercooling (superheating) denotes a process with 
increasing (decreasing) h. The starting configuration of 
the undercooling is a collapsed state, while that of the 
superheating is a smooth state in each surface. 1 x lO'' 
MCS were done at every value of 6, and the final config- 
uration obtained at previous h was assumed as the initial 
configuration of the next h in the processes. The obtained 
hysteresis is consistent to the first-order collapsing transi- 
tion, which was confirmed from the double peak structure 
in h{X'^) in Figs.Etg) and[6th). 

Now we turn to the transition of surface fluctuations. 
The bending energy S^/Nb is plotted in Figs. Wi&) and 
[9jb) against 5, where Nb is the total number of bonds ex- 
cluding the boundary bonds of the holes; S2 is defined only 
on the internal bonds. The variation of S2/NB against b 
becomes rapid with increasing N and hence is slightly dif- 
ferent from that of the model with large-sized holes in [22] , 
where no transition of surface fluctuations was observed. 

The specific heat Cs2 is defined by the variance of the 
bending energy 6*2 such that 



Cs.^^^{iS2-{S2)f), 



(10) 



X' 



N2/H 



(8) 



which is expected to refiect the transition of surface fiuctu- 
ations. Figures [TUfa) andfWb) show Cgj versus b. The ex- 
pected anomalous peaks C'g^^ are apparently seen in the 
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Fig. 9. The bending energy S2/NB versus b of the surfaces (a) 
L = 3and (b) L = 4. 
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Fig. 10. The specific heat Csj versus h of the surfaces (a) 
L = 3 and (b) L = 4. Anomalous peaks indicate a transition of 
surface fluctuations. 



figures, and this indicates the transition of surface fluctua- 
tions because the height of C^^^ increases with increasing 
A'^. In the case of the model with small sized holes, we know 
that C^'^^ decreases with increasing N [22 . Therefore, the 
anomalous structure of Cg^ distinguishes the model in this 
paper and that in [22] , although the size of holes is the only 
difference between the two models. 




5000 10000 
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N 



5000 10000 
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Fig. 11. The log- log plots of Cs^"" against TV of the surfaces 
(a) 1/ = 3 and (b) L = 4. The largest three data in (b) are used 
in the fitting. 



In order to see the dependence of C^!^^ on N more 
convincingly, we plot C™'*'^ versus iV in a log-log scale in 
Figs. [TlTa) and fTUb) . We see that the data satisfy the 
scaling relation 

C'^^^ ~ N^^ (11) 

where cr is a critical exponent. The straight lines drawn 
on the data in Figs. ITTT a) and [TTT b) are obtained by the 
least squares fitting of the data. The largest three data in 
FigfTlTb) were used in the fitting. The results we obtained 



CT = 0.44 ±0.04 (L = 3), 
cr = 0.44 ±0.15 (i = 4). 



(12) 



and these indicate a continuous transition in both of the 
cases L — Z and L = 4 because both of the results obviously 
satisfy cr < 1. We should note that the possibility of the 
first-order transition is not completely eliminated. In fact, 
the surface size seems insufficient even with N — 16062 
because of the finite size effects. 
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Fig. 12. The transition points h^\CQ{N)\ vs. 1/iV for (a) Q = 
X and (b) Q — S2. The straight lines are drawn by fitting the 
data bc[CQ{N)] as a linear function of 1/A'^. 



In the previous section, we mentioned about the finite- 
size effect in the model of this paper. In order to see this 
in more detail, we show the transition points bc[Cx'^{N)\ 
and hc[Cs^{N)] against \/N in Figs.[T2i;a) andff^b). The 
symbol 6c[Cq(A^)] denotes the value of b where the vari- 
ance Cq {N) of the quantity Q has its peak. The straight 
lines are drawn by fitting the data hc[CQ{N)] as a linear 
function of 1/N . Then we obtain the quantities 6c[C'q(oo)] 
in the limit of A^ — > 00 such that 

bc = lAQ (i = 3), 6c = 1-61 (i = 4), (13) 

where both be are common to the quantities Q = X^ and 

Figures fTSTa') andflSFb) show a scaling property of X^ 
with respect to 1 — 6/6c of the surfaces L = 3 and L — 4, 
where be is the transition point in the limit of A^ — > 00 
shown in Eq. (J13p . Crossover behaviors are almost visible in 
the figures, however, we concentrate on the smooth phase, 
which is more clear than the collapsed phase on the figures. 
The straight lines are obtained by fitting the data X'^ in 
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Fig. 13. The log-log plots of X^ against t = l-b/hc of the 
surfaces (a) L = 3 and (b) L = 4, where &c denotes the ones in 
Ea. l|13p . The straight lines in (a) and (b) are drawn by fitting 
some of the data obtained in the smooth phase close to the 
transition point bc{N) of the surface of size A'^; X^ ~ t~". The 
log-log plots of C^J" against t of the surfaces (c) L — 3 and (d) 
L = 4, and the log-log plots of Cgif" against t of the surfaces 
(e) I/ = 3 and (f) L = 4. The straight lines in (c)-(e) were drawn 
by fitting all of the data, and the line in (f) was obtained by 
using the largest three data; Cq'''' ~ i"''«, Q = X^ or Q^S2. 



the smooth phase close to the transition point bc{N) in the 
surface of size N such that X"^ ~ i~", where t= 1 — b/ be- 
Then, we have the exponents a such that aL=3 = 0.76 ± 
0.04 and aL=4 = 0.74±0.09. The two exponents are almost 
identical within the errors. 

The peak values C'^2^ and Cg"^^, shown in Figs. S] and 
[TOl can also show the scaling relation Cq^^ ^ t~^'^. We 
show the log-log plots of C'^2^ against t in Figs.fTWc) and 
[T3td) and those of Cg^/'^ against t in Figs.[T3te) and[T3tf). 
The exponents /3jf 2 of the straight lines in Figs. [TST c) and 
[TSTd) can be compared to ly in Eq. ([7]) and are consistent 
to the first-order collapsing transition, and Ps2 in Figs. 
flSTe) and \W[ i) can be compared to a in Eq. (|12p and 
are also consistent to the continuous transition of surface 
fluctuations. The fitting in Fig. fOTf) was done by using 
the largest three data. Two exponents in Figs. fTSF c) and 



flBT d) are identical within the errors, and those in Figs. 
\Wi e) and [TWf) are also considered to be identical within 
the errors. 

Finally, we comment on the value of Si/N, which is 
expected to be 3(7V— l)/(2A^)~3/2. Our simulation data 
satisfy this relation, and hence the simulations are consid- 
ered to be performed successfully. 



4 Summary and Conclusion 

In this paper, we have studied the conventional surface 
model on triangulated spherical lattices with many fine 
holes, whose size is assumed to be negligible compared to 
the surface size in the thermodynamic limit. The purpose 
of the study is to see whether the phase structure is de- 
pendent on the size of holes or not in the surface model 
with many holes. 

Two types of surfaces are investigated: The first is a 
spherical surface with holes of size characterized by L = 3 
in the unit of bond length, and the hole size corresponds 
to 6 triangles and hence is of hexagonal shape. The sec- 
ond is a surface with holes of size L = 4, and the hole 
size corresponds to 13 triangles. These holes are the min- 
imum size and the second minimum size in the surfaces 
constructed in this paper and in ^^. Therefore, the size 
of holes in the model of this paper is considered to be neg- 
ligible compared to the surface size in the thermodynamic 
limit, while the size of holes in |22] is comparable to the 
surface size in the same limit. 

We find that the model in this paper undergoes a 
discontinuous collapsing transition between the smooth 
phase and the collapsed phase. Moreover, not only the 
smooth phase but also the collapsed phase is considered 
to be physical because the Hausdorff dimensions remain 
in the physical bound, i.e., H <3, in both phases. This re- 
sult is identical to that of the homogeneous model in [2D] . 
These observations are identical to those in the model in 
[22], and therefore, we conclude that the collapsing tran- 
sition is not influenced by the size of holes. 

The second observation in this paper is that the model 
undergoes a continuous transition of surface fluctuations 
at the same transition point of the collapsing transition. 
It was reported in ^20] that the homogeneous model un- 
dergoes a first-order transition of surface fluctuations. In 
the case of large-sized holes, no transition of surface fluc- 
tuations occurs in the model '22]. Thus, our conclusion is 
summarized as follows: The transition of surface fluctua- 
tions is influenced by holes in the spherical surface model, 
and moreover the order of the transition changes depend- 
ing on the size of holes. 

As mentioned in the Introduction, the free boundary 
introduced by holes seems to influence the property of the 
transition of surface fluctuations. Our speculations about 
the dependence of the order of the transition on the size of 
holes are as follows: It is possible that the potential barrier 
between the two degenerate vacuums is dependent on the 
size of holes, where the vacuum states are both smooth 
spherical and different from each other only by the orien- 
tation. One vacuum state transforms to the other by the 
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symmetry transformation n -^ n. The barrier between 
the two vacuums becomes lower and lower with increasing 
size of holes, and the barrier, and consequently the transi- 
tion, disappears when the size becomes comparable with 
the surface size [35]. This seems to be connected to the 
reason why the order of the transition depends on the size 
of holes. 

This work is supported in part by a Grant-in- Aid for 
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